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Abstract 



Contributions from the Breit interaction in atomic-structure calculations 
account for 1.3a of the previously reported 2.5er deviation from the Standard 
Model in the 133 Cs weak charge [S.C. Bennett and C.E. Wieman, Phys. Rev. 
Lett. 82, 2484 (1999)]. The updated corrections for the neutron distribution 
reduce the discrepancy further to l.Ocr. The updated value of the weak charge 
is Qw( 133 Cs) = -72.65(28)cx P t(34) theor . 

The present analysis is a higher-order extension of previous calculation [A. 
Derevianko, E-print physics / 000 1046 1 . 
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Atomic parity-nonconserving (PNC) experiments combined with accurate atomic struc- 
ture calculations provide powerful constraints on "new physics" beyond the Standard Model 
of elementary particles Q. Compared to high-energy experiments or low-energy scattering 
experiments, atomic single-isotope PNC measurements are uniquely sensitive to new isovec- 
tor heavy physics ||. Presently, the PNC effect in atoms has been most precisely measured 
by Wieman and co-workers using 133 Cs @. In 1999, Bennett and Wieman ||] updated the 
value of the Cs weak charge by measuring a supporting quantity, the vector transition po- 
larizability (3, and by re-evaluating the precision of atomic structure calculations |5|U from 
the early 1990s. The determined weak charge Q differed from the prediction J/J of the Stan- 
dard Model by 2.5 standard deviations a. The value of the 133 Cs weak charge from Ref. 0] 
(together with other precision electroweak observables) has been employed in numerous ar- 
ticles. In particular, recent theoretical investigations pL[| interpret this 2.5cr deviation as 
possible evidence for extra neutral vector Z-bosons. 

The main focus of the two previous ab initio relativistic calculations for the atomic struc- 
ture of 133 Cs [|5],^| was the correlation contribution from the residual Coulomb interaction 
(i.e., beyond Dirac-Hartree-Fock level). The purpose of this work is to evaluate rigorously 
contributions from the Breit interaction to PNC in 133 Cs. The previous calculations either 
omitted such contributions |5|], or evaluated them only partially 0]. The present analysis is 
a higher-order extension of my recent calculation ||10|| . It is found that the Breit contribution 
corrects the weak charge by 0.9%, reducing the 2.5a deviation from the Standard Model to 
1.2cr. Including a correction for the neutron density distribution in the 133 Cs nucleus further 
reduces the deviation to 1.0a. Thus the result reported here brings the most accurate atomic 
PNC measurement to date || into substantial agreement with the Standard Model. 



The Breit interaction [] 1| arises due to an exchange of transverse photons between elec- 



trons. Its low- frequency limit, employed here, is given by 

= -— !— (aj • acj + (aci ■ rV,)(<x,- • r^)) 

It is convenient to separate the second-quantized Breit interaction into zero-, one-, and 
two-body parts normally ordered with respect to the core: 

The parity-nonconserving amplitude for the 6S1/2 — > IS1/2 transition in 133 Cs can be 
represented as a sum over intermediate states mPx/2 

_ (7S\D\mP 1/2 )(mP 1/2 \H w \6S) 
^pnc - 2^ p p 

| ^ (7g|^|mP 1/2 )(mP 1/2 |£)|6g) (1) 

rn B 7S - E mPl/2 



Here D [|13| and Hw are electric-dipole and weak interaction matrix elements, and E{ are 
atomic energy levels. It is convenient to break the total Breit correction SEp^c into three 
distinct parts due to corrections in the weak interaction and dipole matrix elements, and 
energy denominators, respectively 

5E PNC = £pnc(5#w) + E PNC {SD) + E P uc(SE) . (2) 

The overwhelming contribution from parity- violating interactions arises from the Hamil- 
tonian 
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G f 

Hw = y|<5wPnuc(^)75 , (3) 

where Gf is the Fermi constant, 75 is the Dirac matrix, and p nnc {r) is the neutron density 
distribution. To be consistent with the previous calculations the p nuc (^) is taken to be a 
proton Fermi distribution employed in Ref. ||. The slight difference between the neutron 
and proton distributions is addressed in the conclusion. The PNC amplitude is expressed 
in units of 10 _11 i|e|ao(— Qw/N), where N = 78 is the number of neutrons in the nucleus of 
133 Cs. In these units the results of past calculations for 133 Cs are -Epnc — —0.905, Ref. 0, 
and -Epnc = —0.908, Ref. ||. The former value includes a partial Breit contribution +0.002, 
and the latter includes none. The reference Coulomb-correlated amplitude 

Ep NC = -0.9075 (4) 

is determined as an average, with the partial Breit contribution removed from the value of 
Ref. [§. 

Hartree-Fock analysis — Before proceeding to the correlated calculations discussed in the 
second part of this work, it is worth examining the Breit contribution to the PNC amplitude 
at the lowest-order level. The conventional Dirac-Hartree-Fock (DHF) equation reads 

Od + Vhf) <Pi = £i<fri , (5) 

where ho is the Dirac Hamiltonian including the interaction of an electron in state i with 
a finite-size nucleus. Vhf is a mean-field Hartree-Fock potential; this potential contains 
direct and exchange Coulomb interactions of electron i with core electrons. A set of DHF 
equations is solved self-consistently for core orbitals; valence wavefunctions and energies 
are determined subsequently by "freezing" the core orbitals. The Breit-Dirac-Hartree-Fock 
(BDHF) approximation constitutes the introduction of the one-body part of the Breit in- 
teraction into the above DHF equation 

(h D + V HF + B«) & = (6) 

Compared to the DHF equations, energies, wave-functions, and the Hartree-Fock potential 
are modified, as designated by tildes. This self-consistent BDHF approximation was used 



by Lindroth et al. [[L2| and a related iterative analysis was considered by Johnson et al. |14 
Both papers point out the importance of the "relaxation" effect, which leads to modification 
of the Hartree-Fock potential through adjustment of core orbitals. In the present work, the 
relaxation effect is taken into account automatically by direct integration of Eq. (Q). 

Most of the Breit contribution to the PNC amplitude can be determined by limiting the 
summation over intermediate states in Eq. ([I]) to the two lowest valence P1/2 states: 6P1/2 
and 7Pi/2- In the DHF approximation one then finds -Epnc = —0.6888 (90% of the total 
value). The lowest-order corrections to matrix elements and energy denominators calculated 
as differences between BDHF and DHF values are listed in Table |. The resultant BDHF 
corrections to -Epnc are: 

E PNC (5H W ) = 0.0022 (0.32%) , 
E PNC (5D) = 0.0020 (0.29%) , (7) 
Epxc(5E) = -0.0019 (-0.28%) . 
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The sum of these three terms leads to dEp^c — 0.0023 in agreement with the 0.002 correction 
found by Blundell et al. ||,[n|. Inclusion of intermediate states beyond 6Pi/ 2 and 7P\/ 2 leads 
to a small additional modification to <5£pnc of -0.00004. Note that if experimental energies 
(which effectively include the Breit interaction) are used in the energy denominators of 
Eq. ([[]), then the £pnc(<5-E) term must be excluded and the total correction becomes twice 
as large: 5E PNC = 0.0042. 

With further examination of the modifications of individual uncorrelated matrix elements 
presented in Table |, one notices the following. 

(i) Weak interaction matrix elements are each reduced in absolute value by 0.3%, which is 
directly reflected in a 0.3% correction to the PNC amplitude. 

(ii) Modification of dipole amplitudes is strongly nonuniform. There are substantial cor- 
rections only to the 6Si/ 2 — 7Py 2 (0.5%) and 7Si/ 2 — 6P1/2 (0.1%) matrix elements. The 
large 0.5% Breit correction to (6Si/ 2 |/}|7Pi/ 2 ) provides partial resolution to a long-standing 
discrepancy of spectroscopic experiment [|nj and ah initio calculations |TT|-|TTf|. The rela- 
tively large Breit correction is caused both by an accidentally small matrix element and by 
admixture into (6Si/2\D\7Pi/2) from a 30 times larger 7Si/ 2 — 7P\/ 2 matrix element. 

(iii) The largest modification in the energy denominators is 0.1% for Ejs — Eqp; however, 
this leads to a 0.3% correction E PN c($E). As recently emphasized by Dzuba et al. pil, such 



large sensitivity of the resulting PNC amplitude to small variations in individual atomic 
properties entering Eq. ([!]) arises due to a cancellation of relatively large terms in the sum 
over states. 

Correlated calculations — It is well known that correlations caused by residual Coulomb 
interactions not included in the Hartree-Fock equations can lead to substantial modifica- 
tions of the lowest-order values. For example, the weak matrix element (6Si/ 2 |-ffw|6-Pi/2) 
is increased by a factor of 1.8 by correlations due to residual Coulomb interactions. It 
will be shown that the correlations are also important for a proper description of the Breit 
corrections. 

The major correlation effects in atoms appear because of shielding of externally applied 
(e.g., electric) fields by core electrons and an additional attraction of the valence electron by 



an induced dipole moment of the core plfl . The former effect is described by contributions 
beginning at second order and the latter in third order of many-body perturbation theory 
(MBPT). Since these two effects lead to the dominant contributions in Coulomb-correlated 



calculations, the third-order analysis reported here seems sufficient |22 



MBPT calculations were performed with the two-body Breit interaction treated 
on equal footing with the residual Coulomb interaction. Sample many-body diagrams are 
presented in Fig. [TJ. To treat the one-body contribution an extension of the B-spline 
basis set technique p3| was developed, based on the Breit-Dirac-Hartree-Fock (BDHF) equa- 
tion (H). Such a formulation made it possible to handle and the associated relaxation 



effect exactly. Contributions of negative-energy states, discussed for example in Ref. |23 



were also included and found to be relatively small |]10| . Two series of third-order calcu- 
lations were performed, first with the Breit and Coulomb interactions fully included using 
the BDHF basis set, and second in the DHF basis set without the Breit interaction and 
negative-energy states. The obtained differences are the Breit corrections reported in Ta- 
ble |. 

Breit corrections to 133 Cs hyperfine-structure magnetic-dipole constants A are discussed 
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first, since these were considered in the literature previously. The correction to hyperfme 
constants is very sensitive to correlations: e.g., Ref. [18| found a numerically insignificant 



modification for Aqs, while Ref. [|19| , |l0f determined the modification to be large (-4.64 MHz), 
and the approach reported here yields +4.89 MHz. In the calculation of Ref. |0| the 
correction was determined as a difference of the BDHF and DHF values, however such 
approach misses two-body Breit corrections of comparable size. In Ref. [[jj!,|l(J a second order 
perturbation analysis was used for the Breit interaction, but the important relaxation effect 
discussed earlier was omitted. The present calculation incorporates all mentioned diagrams 
and is also extended to third order. Using this same calculational scheme, the corrections to 
hyperfme constants for other states of 133 Cs are +1.16 MHz for 7S\/2, -0.51 MHz for 6P1/2, 
and -0.146 MHz for lP\ji- These corrections improve agreement with experiments for the 
ab initio all-order Coulomb-correlated calculations [|18|] to 0.1% for all states except 6P1/2 
where the discrepancy becomes 0.5%. 

Examination of the third-order corrections listed in Table | reveals the significant effect 
of correlations on the Breit contribution. For example, corrections to weak interaction 
matrix elements become three times larger than those in the lowest order. Compared to 
hyperfine-structure constants there is no cancellation of various contributions to the weak 
interaction matrix elements. Using third-order matrix elements and second-order energies 
the following ab initio corrections are determined: Ep NC (5H w ) = 0.0043, Ep NC (6D) = 
0.0035, and E P ^c{^E) = —0.0028. Thus the lowest-order corrections given in Eq. (^) are 
increased. 

To improve the consistency of the calculation, one can combine all-order Coulomb- 
correlated matrix elements and experimental energy denominators tabulated in Ref. || with 
the present third-order Breit corrections. The results are: 

£ PNC (5#w) = 0.0047 (0.5%) , 
E PNC (SD) = 0.0037 (0.4%) . (8) 

The Breit correction in energy-denominators Epnq(SE) was set to zero because the experi- 
mental energies were extensively used in Ref. [§,0. For example, the experimental energies 
were employed in eight out of ten test cases in the scatter analysis of Ref. |J based on 
Eq. (|I|). The total 0.9% Breit correction, 6E PNC = 0.0084, is two times larger than the 
corresponding lowest-order modification, which is rather common in conventional Coulomb- 
correlated calculations. An even larger 2% Breit correction was found in related calculations 
of the electric-dipole-moment enhancement factor in thallium 

Discussion - Combining the calculated 0.9% Breit correction with the reference 
Coulomb-correlated value, Eq.(|]), one obtains the parity- nonconserving amplitude 

£pnc ( 133Cs ) = -0.8991(36) x 10- n z(-Q w /A0 . 

A 0.4% theoretical uncertainty is assigned to the above result following the analysis of 
Ref. M. Since the Breit interaction contributes at the 0.9% level to the total PNC amplitude, 
even a conservative 10% uncertainty in 5-Epnc barely affects the accuracy of -Epno When 
Epnc i s combined with the experimental values of the transition polarizability (3 M and 
EpNc/P 05 one obtains for the weak charge: 

Q w ( 133 Cs) = -72.65(28) cxpt (34) theor . 
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This value differs from the prediction of the Standard Model Q^ 1 = -73.20(13) by 
1.2(j, versus 2.5a of Ref. f3[], where a is calculated by taking experimental and theoretical 
uncertainties in quadrature. This 1.2c deviation is slightly reduced further by taking into 
account corrections for the neutron nuclear distribution in 133 Cs, estimated but not included 
in the final -Epnc of Ref. ||. Recently Pollock and Welliver |^5| determined the relevant 
modification to be AQ^ 1 = +0.11, which reduces the deviation from the Standard Model 
to 1.0 a. 

The present calculation also provides a large Breit correction to the 6Si/ 2 — IP1/2 
electric-dipole matrix element. Using the ab initio all-order Coulomb-correlated value [IB . 
(6-51/21 PI |7Pi/ 2 ) = 0.279, and adding the 0.7% Breit correction of 0.0019, one finds 
(6S1/2I \D\ \7P\/2) = 0.281 in much better agreement with the 0.284(2) experimental 
value [|1(J. The calculated Breit corrections bring most of the ab initio Coulomb-correlated 
hyperfine-structure constants for 133 Cs fl8| into 0.1% agreement with experimental values. 

To summarize, third-order many-body calculations of the contribution of the Breit in- 
teraction to the 133 Cs parity- nonconserving amplitude -Epnc and relevant atomic properties 
are reported. The difference between the present and the earlier calculations || is due to 
additional inclusion of two-body Breit interaction, correlations, and the consistent use of 
experimental energies. The present analysis is a higher-order extension of my recent calcu- 



lation ||i0| . Since the major correlation effects are included, the present third-order analysis 
seems sufficient. The calculations reveal a 0.9% correction to -Epnc leading to a reduction to 
1.2cr of the recently reported 2.5<r deviation ||] of the 133 Cs weak charge from the Standard 
Model value. If corrections for the neutron distribution in 133 Cs nucleus are included, then 
the agreement between the atomic PNC in 133 Cs and the Standard Model stands at l.Ocr. 
Thus the result reported here brings the most accurate atomic PNC measurement to date || 
into substantial agreement with the Standard Model. 
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ences, Office of Energy Research. Part of the work has been performed at Notre Dame 
University during a visit supported by NSF grant No. PHY-99-70666. Calculations were 
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is thankful to W.R. Johnson and M.S. Safronova for useful discussions and H.R. Sadeghpour 
for suggestions on the manuscript. Help with the manuscript and the stimulating interest 
of R.L. Walsworth is greatly appreciated. 
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TABLES 

TABLE I. Breit corrections to matrix elements and energy denominators in a.u.; 
SX, I = Xbdhf — Xdhf, and SX, I + II + III are the differences in the third order of MBPT. 





6*51/2 — 6-P1/2 


6*51/2 — 7i-\/2 


7<Si/2 — 6-f\/2 


7<S*i/2 — 7Pi/2 


H w , DHF 


0.03159 


0.01891 


0.01656 


0.009913 


SH W , I 


-0.00010 


-0.00006 


-0.00005 


-0.000031 


SH W , I+II+III 


-0.00028 


-0.00016 


-0.00014 


-0.000084 


D, DHF 


2.1546 


0.15176 


1.8017 


4.4944 


SD, I 


0.0001 


0.00073 


0.0019 


-0.0004 


SD, I+II+III 


-0.0004 


0.00077 


0.0020 


-0.0012 


AE, DHF 


-0.041752 


-0.085347 


0.030429 


-0.013166 


SAE, I 


-0.000020 


0.000003 


-0.000030 


-0.000007 


SAE, I+II 


-0.000045 


-0.000023 


-0.000034 


-0.000012 
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FIGURES 
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(a) (b) 

FIG. 1. Sample many-body diagrams included in the calculations. Dashed (solid) hor- 
izontal lines represent the Breit (Coulomb) interaction. All orbitals are obtained in the 
Breit-Dirac-Hartree-Fock approximation. Diagram (a) is one of the contributions in the ran- 



dom-phase approximation, and diagram (b) is one of the Brueckner-orbital contributions [21] 
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